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mapped onto a free spinless fermion model with site-dependent hopping amplitudes 
and on-site potentials via the Jordan- Wigncr transformation. Although the single- 
particle eigenstates of this system cannot be obtained in closed form, it is shown that 
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boundary conditions. This proves a conjecture of Hikihara and Nishino [Hikihara T 
and Nishino T 2011 Phys. Rev. B 83 060414(R)] based on numerical evidence. 
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1. Introduction 

The effect of boundary conditions on tlie thermodynamic properties is usually believed 
to be negligible. It is, however, well known that the boundary conditions make a 
difference in the extrapolation procedure from finite-size results to the thermodynamic 
ones. For instance, a leading finite-size corrections to the ground state energy of a 
system with open boundary conditions (OBC) is in general larger than that of a system 
with periodic boundary conditions (PBC). Furthermore, OBC often leads to boundary 
oscillations such as Friedel oscillations. To suppress the finite-size effect due to the 
presence of open edges, the PBC are more often used than the OBC. However, in 
numerical studies, in particular, in the density matrix renormalization group (DMRG) 
method which is a powerful method to study one- dimensional (ID) strongly correlated 
systems, the OBC are more preferable than the PBC. This is because, although there are 
variants of DMRG that can work with PBC [Il[2], they require additional computational 
resources compared with the conventional scheme with OBC. There is a way to remove 
the boundary effects in the DMRG scheme by turning off the interactions smoothly 
around the edges [H H] . 

Recently, a variant of the scheme for the smooth boundary condition, which is 
called the sine-square deformation (SSD), has been proposed |5]. In this scheme, the 
interaction energy scale at the position x in the chain of length L is modified according 
to the following function: 



where < x < L. In [5j, ID free fermion systems with SSD were examined and it was 
found that the SSD completely suppresses the boundary effect in the ground state, i.e., 
the finite-size correction to the ground state energy is identical to that of the chain with 
PBC. Furthermore, the authors found that local correlation functions such as the bond 
strength are translationally invariant, which indicates that the ground states of PBC 
and OBC with SSD are extremely close to each other. A class of generalizations of SSD 
called sinusoidal deformation has been applied to the free fermion systems as well as 
interacting systems such as the extended Hubbard model [6]. It was again found that 
the leading finite-size correction to the ground-state energy has the same form as that 
of the uniform chain with PBC. Motivated by these observations, Hikihara and Nishino 
have studied the SSD in several ID quantum spin systems at criticality such as the XXZ 
spin-| Heisenberg chain in a field [7]. They examined the ground-state properties such 
as entanglement entropy and correlation functions numerically. Their results strongly 
suggest that the ground state of an open chain with SSD is extremely close to that of a 
periodic chain in which the interactions are translation invariant. This conclusion was 
further supported by examining the wave-function overlap, which is very close to unity. 
In particular, for the XY case (A = 0), the wave-function overlap is exactly 1 within the 
numerical accuracy. Note that the deviation of the overlap from the unity becomes of 
the order of 10~^ when A 7^ and hence it is unlikely that the correspondence between 




(1) 
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the periodic chain and the open chain with SSD is exact except for A = 0. 

The aim of this communication is to prove this correspondence analytically for the 
XY case. Although the model can be mapped onto the free spinless fermion system 
via the Jordan- Wigner transformation, hopping amplitudes and on-site potentials are 
position dependent. To our knowledge, the single-particle eigenstates of this free-fermion 
system have not been obtained in closed form, similarly to the Hofstadter problem [8]. 
This is in contrast to the known inhomogeneous XY models whose single-particle states 
can be obtained analytically [9l [101 HI]. Nevertheless, we shall show that the Slater 
determinant of plane waves in the periodic chain is the exact ground state of the open 
chain with SSD. To this aim, we first consider the XY chain with cosine deformation, 
in which the interaction energy scale at the position x is modified by the function 
cos[^(x — |)]. Then we show that the Slater determinant state is a zero-energy state 
of this model. With a simple trigonometric identity, we find that the cosine-deformed 
Hamiltonian consists of the Hamiltonian for the periodic chain and that for the open 
chain with SSD. Since the Slater determinant state is an exact eigenstate of the periodic 
chain, it is obvious that this state is also an eigenstate of the SSD Hamiltonian. From 
the Perron- Frobenius theorem, we can conclude that this state is the unique ground 
state of the SSD Hamiltonian. 



2. Definition of the XY spin chain in a field with sine-square deformation 

The XY spin chain in a magnetic field with PBC is described by the following 
Hamiltonian: 

L-1 L 

where —J < 0, L is the total length of the chain, and 3°' (a = x,y,z) are spin-| 
operators defined on the jth site. This model can be solved exactly using the Jordan- 
Wigner transformation fi2[ [T3] . It has received considerable attention for a long time 
and various quantities such as the correlation functions [HI [151 [HI [IZjj the emptiness 
formation probability [HI [19], and the entanglement entropy [20l [2T1 [22| [23] have been 
calculated exactly. The Hamiltonian with positive exchange J can be obtained from 
"HxY by the unitary transformation U = Ilj odd ^■^P(^^'^i)• 
L-l L 
UHxyU^ = J + S]S]^,) + (-1)^7(^2^1" + SlSf) -hY, S^, (3) 

i=i i=i 
where the case of even L corresponds to the PBC while that of odd L corresponds to 
the antiperiodic boundary condition (APBC) (see footnote 18 of Ref. [7]). 

Let us next introduce the Hamiltonian for the XY chain with the SSD [7|: 

L-1 L 
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where is introduced in Eq. ([T]). As is obvious, this model is defined on the open 
chain of length L, i.e., there is no interaction between the end sites (1 and L). This 
model can also be mapped onto the free spinless fermion system via the Jordan- Wigner 
transformation. However, in contrast to the case of ([2]), the hopping amplitudes and on- 
site potentials result in position dependent. Therefore, the standard Fourier transform 
cannot be applied to solve the problem. 

Before going into details, we give a brief synopsis on the known results of the 
uniform XY chain with PBC, which will be of use later on. The Hamiltonian "Hxy can 
be written in terms of spinless fermions Cj clS 

^XY = ^ J2(^h+^ + sVi^i) + ^r(c[ci + c{cl) -hJ2 - \) ' (5) 
j=i j=i V / 

where the Jordan- Wigner fermions are defined through 

i-i i-i . 

St = c]l[{l-2clcd, S- = c,l[{l-2clc,), S^ = c^c,--, (6) 

i=l 1=1 

with Sf = Sj ± iS'J. The operator T in Eq. (jS]) is given by F := 11^=1 (~'^|) 5 which 
commutes with the Hamiltonian "Hxy- Here, (t| is the Pauli matrix and related to 
Sj through cr| = 25*1. Since F^ = 1, the eigenstates of Tixv are separated into two 
disconnected sectors with F = ±1, in which +/— characterizes configurations with an 
even/odd number of up spins. Therefore, we have 

1 + T 1 — F _ 

"HxY = — ^XY n ^XY5 (''') 



where 



^XY = ~ + ± ^(4ci + c\cl) -hY, Uc, - . (8) 

Now we have to take care of the boundary conditions on the fermions. For F = +1, we 
have to impose APBC on the fermions and they have "half-integer" quasimomenta: 

P=?f^+iV (^eZ), (9) 



L \ 2^ 

where -f < ^ < f - 1 for even L and ^ < ^ < f - | for odd L. We call this 

sector the NS (Neveu-Schwarz) sector in which there is an even number of fermions. On 
the other hand, for F = —1, we require PBC on the fermions and they have "integer" 
quasimomenta: 

p = ^£, (£ez), (10) 

where < £ < ■§ - 1 for even L and -f + ^ < ^ < ■§ - | for odd L. We call this 
sector the R (Ramond) sector in which there is an odd number of fermions. 
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3. Ground state wavefunction of the uniform XY chain 

Our purpose is to show that the ground state of "Hxy is identical to that of "HsXY- We 
shall first summarize the ground state wavefunction of the uniform XY chain. Since the 
magnetization m = '^j=i conserved (in both "Hxy and T-Lsxy), it is convenient 
to work at fixed m. In finite systems, m and the number of fermions (up spins) are 
related to each other through m = ^ — |. In the ground state of TixY, one can show 
that, for a given m, the magnetic field should satisfy 

— Jcos ^^^^ < h < — Jcos ^ — ' ^^^^ 

This gives the equilibrium relation between magnetization and field in the 
thermodynamic limit: 

hm = Jsin(m7r). (12) 

In the XY chain with SSD, we restrict h to the above hm instead of the generic one in 
Eq. ( ITT]) to ensure the solvability of the model. The difference becomes negligible in 
the thermodynamic (L — )■ oo) limit. Note that a similar strategy has been used in the 
numerical studies (see footnote 4 of Ref. [7]). 

For fixed magnetization, the unnormalized ground-state wavefunction of T-Lxv is 
written as 

l^o) = Yl det[e'f'^^]ij=i,2,...,7vki,2:2,...,a:7v) , (13) 

{xi,...,a;jv} 

where \xi,X2, ■■■,Xn) (1 < Xi < X2 < ■ ■ ■ < Xa? < L) denote the configurations of up 
spins at {xi,X2, ...,xn)- Here, the sum extends over all possible spin configurations. In 
the ground state, the set of quasimomenta corresponding to m (A^) is given by 

27r fNl\ 

with i = 1,2,...,N for both the NS and R sectors. Then introducing a new set of 
variables 

Zj := exp ' i^^) 

l^'o) in Eq- ffT3l) is rewritten as 



iV-l 
2 



1*0^= (n^M n {Zi- Zj)\xi,X2,...,XN) , (16) 

{xu...,xn} \i=l / l<i<j<N 

where we have used the property of the Vandermonde determinant. A straightforward 
calculation gives the ground-state energy as 

hL , /L \ sin(7rA^/L) , , 
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Using the magnetization m, this ground state energy can be recast as 



En = —hmL- 



J cosimir) 



hm 



J cosirmr] 



Ti J cosimn) 



+0 



sin(7r/L) \^ ' vr J 6L ' ' 

The leading finite-size correction to the ground-state energy is in agreement with the 
general prediction of the c = 1 conformal field theory [2U [25] . 



4. The XY spin chain with cosine-deformation 

To prove the conjecture of Hikihara and Nishino [7] that the state |^o) defined in Eq. 
(fT6|) is the ground state of "HsXY, it is more convenient to consider the XY model with 
cosine deformation. The Hamiltonian of this system is constructed from Hxy and "HsXY 
as 



T ^ 



cos 



2tt 



E 



cos 



2ti 

T 



We have another decomposition of "He 



Tic 



2 ("^chiral 



chiral / ' 



(19) 



(20) 



where 'Hl^r ■ , is defined by 



± 

chiral 



J 



q-q+ \ „±(i-l/2) qz 



(21) 



with q = exp(i^). We shall consider the eigenvalue problem T-L^ij-^i l"^) = l^') 
instead of that of "Hcos- Note that "H^j^a^i are non-Hermitian and their eigenvalues may 
not be real. For fixed magnetization m and hence fixed number of up spins (A^), the 
eigenstate can be written in the form 

I*) = 5Z '^{Zi,Z2,...,Zn)\Xi,X2,...,Xn) , (22) 

{xi,...,Xff} 

where the coordinates {xi, X2, ...,xn} should be distinct due to the hard-core constraint. 
Let us now derive the Hamiltonians in the first quantized form. To do so, we introduce 
the g-shift operator Tgj which acts on a function of {zj}^^^ as 

). (23) 



Note that F({zj}f^^) is restricted to be a function such that F({zj} 



two out of coordinates coincide. Using Tgj, the action of "H^jj,^! on ^(zi,Z2, ■■■,zn 



is written as i^chirai*(2^i' 2^2, zn) = E^'^^zi, Z2, z^) with 

N , N 



H. 



hr 



chiral 



-1/2, 



N 



N 



H. 



chiral 



(24) 
(25) 
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where we have used the fact that ^^^^ = ^^"=1 <? = 0. Note that the magnetic field 



has been specified according to the relation defined in Eq. f[T2|) . The model is no longer 
solvable even at the field with a small deviation from hm- This is because the Zeeman 
term with the cosine deformation (or its chiral counterpart) does not commute with the 
exchange term with the same deformation. 

We shall show that the Slater determinant state 



JV-l 
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\i=l / l<i<j<Af 

is the exact zero-energy state of i^^jrai- Using a similarity transformation, we have 

[v[/o({.,}f=,)]-^i/,tiral*o({^,}f=l) 

7 ^ / -1 
J \^ / _iV^ -TT — Zk „i TV Q Zj — Zk 

J ^ 



j=l ^ fc(^i) ^ k{^j) ^ 

N 

. f$:(g-/^ + g"-/^)g"^/^.„ (27) 
i=i 

where we have used the relation (fT2l) . Note that \l/o({2j}^i) is nonvanishing for any 
configuration {^jjj^^, which will be shown in the next section. Using the following 
identity 

N . ^ 

E%n^=E^.. (28) 



i^cos^0(^l, 2:2, Zn) = 7T(i^chiral + ^chiral) ^0 (^1 , ^2, ^iv) = 0, (29) 



where t is an arbitrary complex number, we can further simplify Eq. f l27j) . As a result, 
we obtain [^o{{zj})]-'H+^^^^o{{zj}) = 0, which means H+.^^^^i/o{{zj}) = 0. Along the 
same lines as above, we can also show that -f^diii.ai^o({%}) = using the identity Eq. 
(p8l) . Therefore, we arrive at 

1 

2' 

where i^cos is the first quantized Hamiltonian for "Hcos- Therefore, l^'o) in Eq. ( IT6|) is 
the exact zero-energy state of Hcos- It should be noted that |\E'o) is not necessarily the 
ground state of T-Lcos- Since |\E'o) is the eigenstate of TixY with the eigenenergy Eq, we 
immediately arrive at 

^sXy|^o) = yI^o)' (30) 

where we have used the relation among "Hcos, 'Hxy, and "HsXY (Eq- dH]))- Therefore, 
l^'o) is an exact eigenstate of the XY model in the field with SSD. 

Several comments are in order. The remarkable identity Eq. fl25]) can be directly 
proved by comparing the residues of both sides. Another way of obtaining Eq. (128|) is to 
make use of the generating function, which is used to construct Macdonald's g-difference 
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operators whose eigenfunctions are the Macdonald polynomials (see Sec VI in [25]). The 
generating function is defined through 

j=l ■> r>0 

where go{zi, Zn] 0,t) = 1 and 

(7,(^i,...,^„-0,t) = (l-t)V4 n 5^^^ (r>l). (32) 

^-^ ■' Zj — Zh 

j=i k{^j) ^ 

By comparing the coefficients of of both sides in Eq. f l3ip . one can obtain Eq. (12S]) . 
An interesting corollary from the fact K^iirai l^o) = is that |\l/o) is also a zero-energy 
state of the XY Hamiltonian with sine-deformation defined by 

"^sin ~ ^("^chiral "^chiral) 

i=i ^ ^ i=i 
When h = 0, the free-fermion system corresponding to "Hsin is related to the Hofstadter 
(or Harper) problem. This indicates that one can explicitly write down the many-body 
eigenstate of the Hofstadter problem on a thin torus [27] at half-filling even though 
one cannot obtain the single-particle eigenstates in closed form except when the single- 
particle energy is zero [28]. Many other variants can be constructed from 'H^j^.g^i and 
^ciiirai' which |\l/o) IS the exact zero-energy state. It is worth mentioning that the 
explicit definition of q is not necessary for the proof of the zero-energy state. Instead, we 
have used the relation between Zj and Xj, i.e., Zj = q^^ (see Eq. (flSll ). We can, therefore, 
apply the above argument to a large family of models including the systems with the 
hyperbolic deformation [291 EOl EI] and those with the exponential deformation [32] , in 
which the parameter q is taken to be real. 

5. Uniqueness of the ground state 

So far, we have shown that Slater determinant state |\I^o) is the exact eigenstate of 
T-LsXY- In this section, we shall prove that |\l/o) is the unique ground state of "HsXY- For 
fixed magnetization, we specify the basis state we work with as \xi,X2, ■■■,xn), where 
1 < Xi < X2 < ... <a;Ar < L denote the positions of up spins. In this basis, all of the 
off-diagonal matrix elements of %sXY (Eq- (HI)) are non-positive. Since > for all 
j = 1, 2, L — 1, it is obvious that "HsXY satisfies the connectivity condition. Therefore, 
the Perron-Frobenius theorem applies and hence the ground state is nondegenerate and 
can be written as a linear combination of all the basis states: 

(34) 

{xi,...,Xff} 



27r 

TV- 2 



S'. (33) 
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where the coefficients Cx^,x2,-,xi^ are strictly positive for any (xi, X2, xat). Let us 
examine whether the state |\l/o) in Eq. f|T6l) has the same property. Using the relation 
(fT5|) . one can write down the component of |\l/o) as 

^o(^i,--,27v) = (-2i)^^^^ JJ sin{^{xj-Xi)Y (35) 

l<i<j<N 

Therefore, the state 

l^^o) = (-2i)-^^ 1*0) (36) 

can be expressed in the form of Eq. (IMI) . Here, we have used the fact 1 < Xj — Xi < L 
when i < j. Since there can be no other state with positive coefficients only that is 
orthogonal to the ground state, |\l/o) is identical to 1^1/ gs) apart from an overall factor. 
This proves that |\l/o) is the unique ground state of T-LsXY- 

6. Concluding remarks 

In this communication, we have studied the XY chain in a field with SSD. We have shown 
that the unique ground state of this model is identical to that of the uniform XY chain 
with PBC in each sector of fixed magnetization. This explains the previous observations 
in numerical studies that the SSD completely suppresses the effect of open boundaries 
on the ground state of the free-fermion or XY chain [5l [6l [7] . We also found that the 
simultaneous ground state of the uniform and SSD XY chains is the zero-energy state 
of a large class of Hamiltonians including the XY chain with sine/cosine deformation. 
The identity ( l28l) playing a crucial role in the proof has its origin in Macdonald's q- 
difference operators. We can generate an infinite number of identities from Eq. ( 13T]) . 
Therefore, it is interesting to ask whether these identities give a new class of models 
whose ground states can be obtained exactly. In contrast to the XY (free-fermion) 
case, previous works [7] suggest that it is unlikely that the correspondence between 
the periodic chain and the open chain with SSD is exact in the general XXZ chain. 
However, this does not exclude the possibility of another deformed Hamiltonian whose 
ground state is exactly identical to that of the periodic XXZ chain. In this respect, it is 
worth recalling that the ground state of the XXZ chain in the antiferromagnetic regime 
is also an eigenstate of the corner transfer matrix for the six- vertex model [33]. Although 
the corner transfer matrix is well-defined only in the massive regime, the models with 
hyperbolic deformations have some similarities with the corner Hamiltonian [29] . 
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